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Abstract 



^ ; We present measurements of mechanical vibrations induced by 0.6 GeV elec- 
g ! trons impinging on cylindrical and spherical aluminium resonators. To monitor the 
^ " amplitude of the resonator's vibrational modes we used piezoelectric ceramic sen- 
3 sors, calibrated by standard accelerometers. Calculations using the thermo-acoustic 
^ ■ conversion model, agree well with the experimental data, as demonstrated by the 
O ' specific variation of the excitation strengths with the absorbed energy, and with the 
traversing particles' track positions. For the first longitudinal mode of the cylindri- 
^ ^ cal resonator we measured a conversion factor of 7.4 ±1.4 nm/J, confirming the 
^ ; model value of 10 nm/J. Also, for the spherical resonator, we found the model values 
for the L=2 and L=l mode amplitudes to be consistent with our measurement. We 
thus have confirmed the applicability of the model, and we note that calculations 
based on the model have shown that next generation resonant mass gravitational 
wave detectors can only be expected to reach their intended ultra high sensitivity if 
they will be shielded by an appreciable amount of rock, where a veto detector can 
reduce the background of remaining impinging cosmic rays effectively. 
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1 Introduction 



A key issue for a Resonant Mass Gravitational Wave Detector of improved sen- 
sitivity with respect to the existing detectors, is the background due to impinging 
cosmic ray particles []2|, The energy deposited in the detector's mass along a 
particle's track may excite the very vibrational modes that are to signal the pass- 
ing of a gravitational wave. Computer simulations of such effects are based on the 
thermo-acoustic conversion model and earlier measurements of resonant effects in 
Beron et al. and Grassi Strini et al. [j^. According to the model, the energy de- 
posited by a traversing particle heats the material locally around the particle track, 
which leads to mechanical tension and thereby excites acoustic vibrational modes 
At a strain sensitivity of the order of 10^^^ envisaged for a next generation 
gravitational-wave detector, computer simulations |7| show that operation of the 
instrument at the surface of the earth would be prohibited by the effect of the cos- 
mic ray background. Since the applicability of the thermo-acoustic conversion model 
would thus yield an important constraint on the operating conditions of resonant 
mass gravitational wave detectors, Grassi-Strini, Strini and Tagliaferri measured 
the mechanical vibrations in a bar resonator bombarded by 0.02 GeV protons and 
g*]^Q-4 QgY electrons. We extended that experiment by measuring the excitation 
patterns in more detail for a bar and a sphere excited by 0.6 GeV electrons. Even 
though we cannot think of a reason why the model, if applicable to the bar, would 
not hold for a sphere, we did turn to measuring with a sphere also. 
We exposed two aluminium 50ST alloy cylindrical bars and an aluminium alloy 
sphere, each equipped with piezoelectric ceramic sensors, to a beam of « 0.6 GeV 
electrons used in single bunch mode with a pulse width of up to « 2/is, and ad- 
justable intensity of 10^ to lO^'^ electrons. We recorded the signals from the piezo 
sensors, and Fourier analysed their time series. Before and after the beam run we 
calibrated the sensor response of one of the bars for its first longitudinal vibrational 
mode at « 13 kHz to calibrated accelerometers. 

2 Experiment setup and method 

In the experiment we used three different setups in various runs, as summarised in 
table I: two bars and a sphere. With the un-calibrated bar BU we explored the 
feasibility of the measurement. Also, bar BU proved useful to indirectly determine 
the relative excitation amplitudes of higher longitudinal vibrational modes, see sec. 
[4.i| . With bar BC calibrated at its first longitudinal vibrational mode, we measured 
directly its excitation amplitude in the beam. Finally, with the sphere we further 
explored the applicability of the model. 



Table I. Characteristics of our setup. 



Setup code name: 


BC 


BU 


su 


Resonator type: 


bar 


bar 


sphere 


Diameter: 


0.035 m 


0.035 m 


0.150 m 


Length: 


0.2 m 


0.2 m 


- 


Suspension: 


plastic string 


plastic string 


brass rod 0.15 m*0.002 m 


Piezo sensors: 


1 


2 


2 


Piezo hammer: 





1 


1 


Capacitor driver 


1 








Direct calibration 


yes 


no 


no 


Beam energy 


0.76 GeV 


0.62 GeV 


0.35 GeV 


Beam peak current 


3 mA 


18 mA 


19 mA 


Electrons per burst 


« 10'' 


^ 5 ■ 10"' 


^ 5 ■ lO^" 


Mean absorbed energy per electron 


0.02 GeV 


0.02 GeV 


0.1 GeV 


Typical absorbed energy per burst 


0.01 J 


0.6 J 


3.0 J 



2.1 Electron beam 

We used the Amsterdam hnear electron accelerator MEA [TDj delivering an elec- 
tron beam with a pulse-width of up to 2 fis in its hand-triggered, single bunch 
mode. The amount of charge per beam pulse was varied, recorded by a calibrated 
digital oscilloscope, photographed and analysed off line to determine the number of 
impinging electrons per burst. 

2.2 Suspension and positioning 

In both setups BC and BU, see fig. 1, the cylindrical aluminium bar was horizontally 
suspended in the middle, as indicated in the figure, with a plastic string. The bar's 




Figure 1: Calibrated bar setup BC. The beam traverses the bar's front side perpendicular to the drawing 
plane. 

cylinder axis was positioned at 90^ to the beam direction. The bar's suspension 
string was connected to a horizontally movable gliding construction, enabling us to 



handle the resonator by remote control, and let the impinging electron beam hit it at 
different horizontal positions. The aluminium sphere SU, see fig. 2, was suspended 
from its centre by a brass rod. Either the bar's gliding construction, or the sphere's 
suspension bar, was attached to an aluminium tripod mounted inside a vacuum 
chamber [)^, which was evacuated to about 10^^ mbar. By remote control, we 
rotated the tripod and moved it vertically to either let the beam pass the resonator 
completely, or let it traverse the resonator. We let the beam traverse the sphere 
at different heights and different incident angles with respect to the piezo sensors 
positions on the sphere. We mark the beam heights as E (Equator) and A (Africa) 



suspension rod 




Figure 2: The spherical resonator setup, SU. 



at 0.022 m below the equator. The E beam passed horizontally through the sphere's 
origin, remaining in the same vertical plane for the A beam. 

2.3 Sensors and signal processing 

In setup BC we used a single piezo sensor of ~ 15 * 3 * 1 mm^ and glued it over it 
full length at 0.01 m off the centre on the top of the bar. Bar BC was equipped with 
a capacitor plate of 0.03 m diameter at a distance of ^ 0.004 m from one of its end 
faces. In setup BU one piezo sensor of ~ 3 * 6 * 0.3 mm^ was fixed on one end face 
of the bar. 

A similar sensor of about the same dimensions was fixed in the same manner, ori- 
ented parallel to the cylinder's long axis at a position 35 mm away from the end-face. 
In the third setup, SU, see fig. 2, two piezo sensors of ^ 3 * 6 * 0.3 mm^ were glued 
to the sphere's surface. One was situated at the equator, with respect to the vertical 
rotation axis, the other one at a relative displacement of 45^ west longitude, and at 
45^ north latitude. 

For the setup in use, each sensor was connected to a charge amplifier of « 2*10^^ V/C 
gain. The signals were sent through a Krohn-Hite 3202R low-pass 100 kHz pre-filter. 



to a R9211C Advantest spectrum analyser with internal 2 MHz pre-sampling and 
125 kHz digital low-pass filtering. The oscillation signals were recorded for 64 ms 
periods at a 4 /is sample rate. The beam pulse could be used as a delayed trigger to 
the Advantest. Using the memory option of the Advantest, the piezo signals were 
recorded from 0.3 ms onward before the arrival of the trigger. The data were stored 
on disk and were Fourier analysed off line. 

2.4 Checks and stability 

The data were taken at an ambient temperature of ~ 23 ^C. By exciting the res- 
onator with the piezo-hammer we checked roughly its overall performance. As to 
be discussed in section 0, setup BC was calibrated before and after the beam run. 
The instrument's stability was checked several times during the run by an electric 
driving signal on its capacitor endplate. 

3 Calibration of bar BC's piezo ceramic sensor 

A standard accelerometer mounted on the bar damped the vibrations too strongly 
to confidently measure their excitations in the electron beam. Therefore the re- 
sponse of the piezoelectric ceramic together with its amplifier was first calibrated 
against two 2.4 gramme Bruel&Kjaer 4375 accelerometers glued, one at a time, to 
bar BC's end face and connected to a 2635 charge amplifier. The resonator was 
excited through air by a nearby positioned loud-speaker driven from the Advantest 
digitally tunable sine-wave generator. The output signals from both the piezoelec- 
tric ceramic amplifier and the accelerometer amplifier were fed into the Advantest. 
Stored time series were read out by an Apple Mac 8100 AV, running Lab- View for 
on-line Fourier analysis, peak selection, amplitude and decay time determination. 
We took nine calibration runs varying the charge amplifier's sensitivity setting, and 
dismounting and remounting either of the two accelerometers to the bar. For the 
lowest longitudinal vibrational mode we calculated the ratio of the Fourier peak 
signal amplitudes, from the piezoelectric ceramic and accelerometer. 
With the calibrated bar BC positioned in the electron beam line we checked the 
stability of the piezoelectric ceramic's response intermittently with the beam runs 
by exciting the bar through its capacitor plate at one end face, electrically driving it 
at and around half the bar's resonance frequency. We found the response to remain 
stable within a few percent. 

After the beam runs we took additional calibration values in air with a newly ac- 
quired Bruel&Kjaer 0.5 g 4374S subminiature accelerometer and a Nexus 2692 A0S4 
charge amplifier. In fig. 3, typical frequency responses are shown when driving the 
bar by a loud-speaker signal. The upper part gives the Fourier peak amplitude 
of the bar's 13 kHz resonance as measured with the accelerometer. The lower part 
gives the corresponding amplitude for the signal from the piezoelectric ceramic. The 
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Figure 3: Response at the lowest longitudinal acous- figure 4: Response of bar EC's piezo (lower) and ac- 
tio frequency of bar EC's piezo (lower) and accelerom- celerometer (upper) Fourier amplitudes by electrostatic 
etcr (upper) Fourier amplitudes by constant amplitude capacitor plate driving. Left: on-resonance, /=6481.7 
loud-speaker driving. Left: on-resonance, /=12950 Hz. ^z- Right: slightly off-resonance, /=6480.0 Hz. The 
Right: slightly off-resonance, /=12850 Hz. Pi^^" P^^'^k of constant amplitude at 6.5 kHz arising from 

crosstalk is absent in the accelerometer, while the acous- 
tic resonance is clearly seen at 13 kHz in both. 

right hand side of the picture shows the ampUtudes to be smaller, as expected when 
driving the bar slightly off resonance. We calculate the decay time, r, of the k-th 
mode amplitude Ak{t) = Ak{0) • e"*/"^ to be r = 0.4 s for this setup, equipped with 
the relatively light accelerometer. 

Figure 4 shows the corresponding two signals when driving the bar by the capacitor 
plate at 6.5 kHz, that is at half the bar's resonance frequency. Here, the direct elec- 
tric response of the piezoelectric ceramic's signal to the driving sine- wave is present, 
clearly without a mechanic signal, as would have shown up in the accelerometer. 
The direct signal at 6.5 kHz remains constant. On the other hand, the bar's me- 
chanical signals on and off its resonance frequency around 13 kHz show the expected 
amplitude change again, thereby demonstrating that around the bar's resonance, the 
piezoelectric ceramic does only respond to the mechanical signal, not to the electric 
driving signal. See also the caption of fig. 4. 

We calculated the average value of Rq = V^^"^'^ /V^^"^'^ and the error over all 
29 measurements, finding for the calibration factor at f=13 kHz, 

^ = RoS{27rff = (2.2 ± 0.3) V/nm (1) 

where S = 0.1 V/ms~^ is the amplifier setting of the accelerometer. 



4 Beam experiments 



Sensor signals way above the noise level were observed for every beam pulse hitting 
the sphere or the bar. We ascertained that: a) the signals arose from mechanical 
vibrations in the resonator, and b) they were directly initiated by the effect of the 
beam on the resonator, and not arising from an indirect effect of the beam on the 
piezo sensors. Our assertion is based on a combination of test results observed for 
both the bars and the sphere, as now to be discussed. 
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Figure 5: The 2 /xs electron bunch (middle trace) and Figure 6: A typical Fourier spectrum of bar BC as 
the piezo sensors signals on the sphere (upper and lower excited by the electron beam. Data were analysed over 
traces). 0.016 s duration, from 0.008 s onwards after the beam 

passed. Identified vibrational modes are indicated. 

First, when the beam passed underneath the resonator without hitting it, we 
observed no sensor signal above the noise. Second, as shown in fig. 5, the sensors' 
delayed responses after the impact of the beam agreed with the sound velocity. Here 
the beam was hitting the sphere at a position 5 mm above the sphere's south pole. 
The middle trace shows the beam pulse of 2//s duration. The two other traces 
show both piezo sensors to respond with a transient signal right from the start time 
of the beam's arrival and to begin oscillating after some delay, depending on their 
distance from the beam. The distance of the equatorial sensor to the beam hitting 
the sphere at the south pole was 0.11 m, corresponding to « 22 ^us travel time for 
a sound velocity of 5 * 10^ m/s. The signal is indeed seen in the lowest trace 
starting to oscillate at that delay time. The upper trace shows the signal from the 
second sensor situated on the northern hemisphere at 0.14 m from the traversing 
beam, correspondingly starting to oscillate with a delay of 28 //s after the impact 
of the beam. Third, after dismounting the piezo-hammer from the resonator, we 
observed that the sensor signals did not change, which showed that the activation 
is not caused by the beam inducing a triggering of the piezo-hammer. Fourth, to 



simulate the electric effect of the beam pulse on the sensors, we coupled a direct 
current of 60 mA and 2.5 fis duration from a wave packet generator to the bar. 
Apart from the direct response of the piezo-sensor during the input driving wave, 
no oscillatory signal was detected above the noise level. Finally, we measured the 
dependence of the amplitudes in several vibrational modes on the hit position of the 
beam, as will be described in the following sections. We found the amplitudes to 
follow the patterns as calculated with the thermo-acoustic conversion model. 



4.1 Results for the bar 



In fig. 6 a typical Fourier spectrum of bar BC is shown up to 55 kHz. The arrows 
point to identified vibrational modes [TTl . From a fit of ii' and /o of the longitudinal 
frequencies fi = L • /o • (1 — L'^K) [jT^ of the modes for L=l,..,4, we find /o = 
12933, K = 0.0022, where /o is related to the sound velocity by Vg = 21 * fo = 
5173 m/s for our bar length of /=0.2 m. For the Poisson-ratio a = 2l^JK / [nr)^ r 
being the cylinder radius of the bar, from our fit we get a = 0.338. The values agree 
well with the handbook [|T3j quoting cr=0.33 and i;s=5000 m/s for aluminium. The 
root mean square error of the fit is 35 Hz, in correspondence with the 30 Hz frequency 
resolution used in the Fourier analysis. Other peaks correspond to torsional and 



transverse modes [jTT], [T^ 
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The Fourier amplitudes Ak of f{t) = E^/ce*"^*^^ of the modes depend linearly (as 
shown for the 13 kHz, L=l mode in fig. 7) on the integrated charge in the beam 
pulse for a fixed beam position, and therefore also linearly on the energy deposited 
by the beam, which ranged in these runs from 0.06-0.8 J. The spread in the ratios 



of the amplitudes to the beam charge, shows the Fourier amplitudes to reproduce 
within ± 10%. 

The agreement of the model to within 10% with the measured data is shown in 
fig. 8. The figure shows the measured Fourier amphtudes of bar BC at the piezo 
sensor and the calculations following Grassi Strini et al. [§, [Tlj as a function of the 
hit position along the cylinder axis for the first four longitudinal modes. For each 
mode the average model value was scaled to the average measured value. The best 
fit was found with a shift of the hit positions along the bar, by an overall offset 
of a:o=-0.0075 m, which corresponds to the crude way we aligned the bar with the 
beam line. 

4.1.1 Lowest bar mode excitation amplitude 



For the 13 kHz, L=l mode we determine the absolute amplitude for a comparison 
with the model calculation of ref. [0, [Tlj. Firstly, we use the amplitude function 
Bo(x), see eq. 9 of ref. [0, [Tlj, by rewriting it in the form: 

Bq{x) = 2-ko- AE/tt X cos{7rx/l)sin{7rT]/ (21)) /tct]/ {21) (2) 

with 

kq = a ■ l/{cy ■ M) = a/{cy- pO). (3) 

In this expressions x is the hit position along the the cylinder axis, / the bar length, 
T] the beam diameter, a the thermal linear expansion coefficient, p the density, c^, the 
specific heat, O the cylindrical surface area of the bar, and AE" the energy absorbed 
by the bar. From Bq{x) we derive the functional form for the measured values of 

where dE/dQ is the beam energy absorbed by the bar per unit of impinging beam 
charge, (3 the calibration factor as discussed in section 0, and D the decay factor 
e~*/^, since eq. applies at excitation time and we have to correct the amplitude at 
measuring time for the mode's decay, corresponding to its Q-factor. Therefore, 

Wsens{x) = Kexp2/7r X COs{tIX / l)sin{'KT] / {21)) / TIT] / {21) , (5) 

with 

dE 

From fitting eq. to the measured values Wsens{x) given in table II with K^xp as the 
free variable, we find our presently measured value for Kq^^ = Kexp ■ /{dE/dQ- D- (3) 
which we compare to the model value in eq. 0. Secondly, the decay time was 
measured by recording the sensor signals after a trigger delayed by up to 1.6 s at a 



fixed beam hit position. An exponential fit A{t) = * e~*^^ to ttie mode amplitude 
gives r = (0.36 ± 0.01) s for the L=l mode. This corresponds to a Q-value of 
~ 15000, a value consistent with the room temperature measurement of aluminium 
as in ref. []T5| , and indicating a negligible influence of the suspension and piezoelectric 
ceramic sensor for this mode. From the measured value of r and a mean delay time 
from the start of the beam pulse of 0.016 s, we calculate the decay factor to be 
D = 0.95. 

Table II. Excitation values Wsens, equalling the ratio of the measured Fourier 
amplitude and the measured beam pulse charge at each of the indicated hit 
positions on the bar for the 13 kHz, L=l mode. 



hit position x cm 





1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


W«en. V/nC 


0.185 


0.216 


0.167 


0.180 


0.225 


0.152 


0.152 


0.157 


0.112 


0.089 


0.057 



Thirdly, as indicated, we use the data for Wsens in the second row of table II to 
fit the variable K^xp in eq. |5|, where now x is the hit position as given in row 1, 
1=0.2 m , and r/ = 0.002 m. The value found in the fit is K^xp = (0.300 ± 0.025) 
V/nC. Fourthly, from a Monte Carlo simulation at the beam energy of 570 MeV 
used for these runs, we calculate the mean absorbed energy and the mean energy 
spread which results from the fluctuating energy losses of the passing electrons and 
the energies of the secondaries escaping from the bar, as AEe=(19± 2) MeV. The 
electron beam pulse thus deposits dE/dQ = (0.019 ± 0.002) J/nC in the bar. Using 
the measured cahbration value at f=12986 Hz as given in eq. [1|, /5=(2.2 ± 0.3) V/nm, 
we arrive at 



4'^'' = (7.4 ± 1.4)nm/J. 



(7) 



Finally, we calculate the model value of «:o from the material constants as being 
HiQ=10 nm/J, neglecting the much smaller error as arising from some uncertainty in 
the parameters. We conclude that Kq^^/kq = (0.74 it 0.14), a result that is consistent 
with the validity of the model of ref. [j^, |T^ . 

The measured maximum excitation amplitude at beam position a:=0, see fig. 8 for 
the 13 kHz, L=l longitudinal mode thus corresponds to (0.13 ± 0.02) nm. 



4.1.2 Higher bar mode excitation amplitudes 

Having determined the correspondence between the model calculation and the ex- 
periment's result for the first longitudinal vibrational mode amplitude, we return to 
some of the higher vibrational modes. To compare the modes we need to take the 
sensor position on the bar into account. We rewrite the displacement amplitude of 
eq. 5 from ref. as a function of hit position Xh and sensor position Xg as: 

^odd-L = {2K/L7T)sin{L7TXs/l)cOs{L7TXh/l), 

^even-L = {2k/ L7T)cOs{L7TXs/l)sin{L7TXh/l) , (8) 



where / is the bar length. We dropped the beam width correction term which 
would lead to a less than 0.1% correction even for L=4. We approximate the sensor 
response by the local strain along bar EC's cylinder axis, that is to the d$/dxs of 
eq. |S|, arriving at a sensor response, S^: 

Sodd-L = BLC0s{L7TXh/l), Bodd-L = {2eti / 1)cos{Lttx s / 1) , 

Seventh = BLSin{L7TXh/l) , B^ven-L = {2eK/l)sin{L7TXs/l) , (9) 

where e is a sensor response parameter. The Xs dependent term did not enter into 
the calculation of Kq^^ in the previous section, since the calibration was done at the 
same sensor position as the beam measurement. However, for a comparison between 
the modes, the dependence on the sensor position Xg has to be taken into account. 
Since the variables are strongly correlated, we, first, fitted for each mode the term 
Bl in the Xh dependent part of eq. |9| to the measured value of Wgens for the mode, 
shifting the origin of Xh by 0.0075 m, as mentioned before. The results are given in 
the first row of table III. 

Table III. Bar BC modes comparison. The piezoelectric ceramic sensor responds 
to the bar's strain. 



description 


symbol 


13 kHz,L=l 


25.6 kHz,L=2 


38 kHz,L=3 


50 kHz,i=4 


amplitude 




0.12±0.01 


0.021±0.002 


0.033±0.03 


0.052±0.005 


decay correction 


D 


1.04±0.001 


1.17±0.02 


1.49 ±0.06 


1.14±0.01 






0.12±0.01 


0.025±0.003 


0.049±0.005 


0.059±0.006 


sensor position factor 


Pl 


1.02 


2.61 


1.20 


1.41 


6^2/; = Bf' X PL{arb.u.) 


t 

K 


0.12±0.01 


0.065±0.007 


0.059±0.006 


0.083±0.008 



Second, we corrected the amplitudes B^^°'^ for the mode decay with a factor 
given in row 2, and corresponding to the times ri = 0.36 s, r2 = 0.10 s, ra = 
0.04 s, r4 = 0.12s, which leads to the values of B^j^^ in row 3. Finally, we multiphed 
with the factor Podd-L = 1/cos{L7tXs/1), Peven-L = 1/ sin{L7TXs/l), where the bar 
length is /=0.2 m. Since the sensor extends from 0.005 through 0.020 m from the 
center of the bar, we use the mean sensor position Xs ~ 0.0125 m. The resulting 
values of = 2€hi/l^ shown in the last row, should be independent of L. For L=2,3,4 
they are rather closely scattered around a mean value of n' = 0.07 which is, however, 
at about half the L=l value. This discrepancy might have originated from some 
resonances of the sensor itself, and we suspect the strong peak at 23 kHz, shown in 
figure 6, to be an indication of such resonances playing a role. 
Since the amplitudes of the higher modes for bar BC do not comply with our ex- 
pectations we turn, as a further check, to our un-calibrated measurements with bar 
BU. It had been equipped with a piezoelectric sensor at one end face where the 
longitudinal modes have maximum amplitude. The sensor had been mounted flatly 
with about half of its surface glued to the bar, and responding to the bar's sur- 
face acceleration, not its strain as at bar BC. We extract the values from our 
measurement analogously as for bar BC, following again the model calculations of 
Grassi Strini et al. []5|, using the L=l mode as the reference. The results are given 
in table IV. 



Table IV. Bar BU modes comparison. The piezoelectric ceramic sensor responds 
to the bar's acceleration. The value of /^^^"^ for the L=l, 13 kHz mode is used 
as the reference for the higher modes. 



description 


symbol 


13 kHz,L=l 


25.6 kHz,L=2 


38 kHz,L=3 


50 kHz,_L=4 


relative amplitude 




1 


1.15 


11.5 


3.6 


relative decay correction 


D 


1 


3±1 


0.7±0.2 


1.3±0.9 


(a;z,=i/a;i)^ 


n 


1 


0.26 


0.12 


0.07 


gmea, x £) X 


meas j meas 


1 


0.8±0.3 


0.9±0.4 


0.3±0.3 



After applying the decay correction factor D and the frequency normalisation 
factor the results should be independent of L. The L=4 value is significantly 
low, which, again, might be due to some interfering resonance. The L=2 and L=3 
values, however, do not significantly deviate from the L=l value, thus confirming 
the model calculations for these higher modes too. 



4.2 Results for the sphere 



Our measurements on the sphere consisted of a) hitting the sphere with the beam 
at one of two heights in the vertically oriented plane through its suspension: at the 
equator (E) and at 0.022 m southward (A); b) rotating the sphere with its two fixed 
sensors over ISO'^ around the suspension axis at each beam height, and measuring 
several times back and forth by steps of 30^ to diminish the influence of temperature 
and beam fluctuations, ending up on a 10*^ angular lattice. The Fourier amplitude 
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Figure 9: The Fourier amplitude spectrum of the Figure 10: Decay of sphere's vibrational modes as mea- 
sphere SU averaged over all measured angles for sensor- sured with delayed data taking of the spectrum analyser. 
1, at beam height position E. The modes and frequencies The curve shows the fit of A{t) = * e^*/^. Upper: 
as calculated, are indicated. 17.6 kHz, L=2, t=1.1 s. Middle: 24 kHz, L=l, r=0.4 

s. Lower: 37 kHz, L=0, r=0.1s. 



spectrum of sensor-1, averaged over the angular positions, is shown in fig. 9. The 
lowest spheroidal mode is most relevant for a spherical resonant mass gravitational 



wave detector, and we therefore focus on a few spheroidal modes. As expected, 
the lowest spheroidal L=2 mode is seen at 17.6 kHz, the lowest spheroidal L=l 
mode at 24 kHz, and the lowest spheroidal L=0 mode at 37 kHz. Some other peaks 
are also indicated in the figure, though not the toroidal modes, which we neglect 
completely. It should be noted that while the L j^O amplitudes oscillate over the 
angles, the L=0 amplitude does not, leading to a relative enhancement of the latter 
in the angle-averaged fig. 9. 

The Fourier amplitudes, again, showed a linear dependence on the deposited energy. 
To determine the decay times at /=17.6, 24 kHz and 37 kHz, see fig. 10, we took 
data with up to 4 s delay in the spectrum analyser, and found r 1 s, 0.4 s and 
0.1 s respectively. 

The angular distributions for the amplitude of the 37 kHz, L=0 mode at the two 
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Figure 11: 37 kHz Fourier- modulus angular distribu- 
tion. The beam hits the sphere at upper: E, lower: A. 

vertical beam positions E and A are shown in figures 11. The L=0 amplitude is 
independent of the angle, and since the amplitude is constant to within 20% we 
infer from the Fourier-modulus' deviation from fiatness a 20% variation of the beam 
intensity from one shot to another. 

During our measurement with the sphere we were unable to use the beam pulse as 
a trigger, implying that the start time of data acquisition with respect to the beam 
pulse is unknown. In our further analysis we will therefore use only the Fourier 
modulus, and not analyse the phases. 

The absolute scale of the 37 kHz Fourier amplitude turned out to be 5 times 
larger than the model value for sensor-2 and « 50 times for sensor- 1. We assume 
this discrepancy to be based on some interference effects, possibly with a sensor 
resonance and a suspension bar mode, and we do not further analyse the L=0 



mode. To disentangle the angular distributions in general, we felt, would squeeze 
the results of our simple measurement too much, for a couple of reasons. First, the 
Fourier amplitude Al for any multi-pole order L in the sphere's case is actually 
a sum of M-submodes. Though they would be degenerate for an ideal sphere, 
in practice some M-modes might or might not turn out to be split beyond the 
frequency-resolution of A/=30 Hz. Second, both sensors si,S2 should be taken to 
have unknown sensitivities, e^^., in three orthogonal directions, with phase factors 
+1 or -1 for their orientation. Third, though each mode would start to be excited 
within the same sub-nanosecond time interval of the beam crossing, the building up 
of each mode's resonance vibration may lead to a specific phase depending on 

the mode's spatial relation to the beam path. 

We now show that the calculated angular distributions have the signature of the L- 
character of the measurement. Therefore, we write the Fourier modulus at different 
impinging beam positions 6^ as a function of the angle (j) as: 

Al,s,h{4>) =1 Fl^s.H ^sj ■ E Si,M,.,6,UL,M^,,^,5,(0)e'^^ I, (10) 

-Ml 

where Fi^sjfik is a frequency response function for each sensor, which may depend 
also on the beam position. This normalisation factor is expected to be of order 
1, and is kept fixed at 1 for the L=2 distributions. It is used as a free parameter 
for the L=l distributions to compensate for the rather inaccurate knowledge of a) 
the sensor positions on the sphere's surface, b) the beam track location and c) the 
electrons and photons shower development along the track, since the exact excitation 
strengths of the modes are quite sensitive to such data. As the first step in the fitting 
procedure we separately calculated the SL^MLA'^L,ML,sj,bk: where sl^MlM mode's 
strength from the beam excitation, as detailed in the appendix. We inserted the 
calculated SL^ML,bk'^L,ML,sj,bk ill ^ hierarchical fitting model, to simultaneously fit []T^ 
the relevant parameters of eq. |rD|to the 17.6 kHz, L=2 Fourier modulus Ai^sj.hk for 
both sensors si and S2 at both beam positions E and A. This fit led to a reduced 
X^=1.3 at 59 degrees of freedom. Next, with fixed values for the sensor efficiencies 
so established, we fitted the relevant parameters for the 24 kHz, L=l Fourier 
peaks, including the L=l sensor response factors F. At all stages the | ^^^,5^ I of the 
phases were kept within the bounds of the period of mode-L. With an uncertainty 
in the beam charge and in the Fourier peak amplitudes of ~ 20% each, the error 
amounts to ~30%, and we took a minimum absolute error of 2x10^^ for sensor si 
and 1x10^^ for sensor S2. In total we have 152 data points, while the total number 
of fitted parameters is 27, including a relative normalising factor for the mean beam 
current at beam position A with respect to the mean current at beam position E. 
We found for the total fit a reduced x^=l.Q at 125 degrees of freedom. The L=l 
response factors remain within 1.1 and 0.2. The results of the fits to the 17.6 kHz 
and 24 kHz are given in figures 12 and 13 and table V. Note the different vertical 
scales used for sensor- 1 and sensor-2 in both pictures. Some of the parameters given 
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Figure 12: Data points (+) and fit results ( — ) for Figure 13: Data points (+) and fit results ( — ) for the 
the sphere's 24 kHz, L=l mode. Left column: sensor-l, sphere's 17.6 kHz, L=2 mode. Left column: sensor-l, 
right colum: sensor-2. Upper row: beam position E, right colum: sensor-2. Upper row: beam position E, 
lower rows: beam position A. The x-axes give the angle lower rows: beam position A. The x-axes give the angle 
of sensor-l. Note that the y-scales are different for the of sensor-l. Note that the y-scales are different for the 
two sensors. two sensors. 



in the table are strongly correlated. 

Table V. Results of a hierarchical fit to the data of the sphere at the 19 mea- 
suring angles of both sensors Si and S2 of the Fourier modulus at 17.6 kHz and 
24 kHz at beam positions E and A. 
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We conclude that the measured Fourier amplitude angular distributions are con- 
sistent with the model value for the L=2 and L=l mode signature. 



4.2.1 The sphere's absolute displacement 

Finally, to estimate the order of magnitude of the sphere's absolute displacement, 
we have to take an intermediate step by first normalising bar BU to the calibrated 
results for bar BC and then use bar BU as a calibration for the sphere. The sensors 
used on bar BU consisting of the same sensor material and having been cut roughly 
to the same size, we assume to be identical to the ones used on sphere SU. The 
amplifiers used are identical. The bar BU sensors, however, differ strongly from 
those of bar BC. 



We arrive at an indirectly calibrated value for eBu,2 = (3.4* 10^"^ ± 30%) V/ms^^ of 
sensor-2 on bar BU. The value of eBu,i for sensor-1 is about ten times smaller. 
Then, for the sphere, the fitted value of eg. given in table V, shows the largest value 
of sensor-2, e^^su = (4 * 10^^ ± 10%) V/ms~2, to lead to a ratio « (1.2 ± 0.4) with 
cbu- Again, the values for sensor-1 are about ten times smaller. The error of ~ 33% 
is the propagated statistical error only. The result seems reasonable. So, the model 
calculation and our sphere measurement results are of the same order of magnitude 
on an absolute scale too. 

From the maximum Fourier modulus, V'^°'^ = (0.003 ± 0.001) V, of the 17.6 kHz, 
L=2 sphere mode measured in sensor sensor-2 as given in fig. 13, and the ab- 
sorbed energy of 3.1 J, we find the maximum sphere's displacement to correspond 
to (0.2±0.1)nm/J. 



5 Discussion 

Having confirmed the thermo-acoustic conversion model in the present experiment, 
we discuss some points about extrapolating these results to the actual operation of 
a resonant gravitational wave detector. Firstly, in our experiment many incident 
particles deposited their energy in the resonator, in contrast to a single muon hit- 
ting an actual detector. However, from this difference it seems unlikely to reach 
different conclusions, especially since in the process of depositing energy along its 
track, the muon will generate lots of secondary particles too. Also, we measured at 
room temperature while actual detectors would have to operate in the millikelvin 
range. An aluminium resonator, for instance, at such a temperature, would be su- 
perconducting, and it is as yet unclear how the decoupling of the electron gas from 
the lattice would affect the process of acoustic excitation. 

Therefore, we hold it of particular importance for the prospected shielding of a next 
generation resonant mass gravitational wave detector that an existing millikelvin 
detector hke the Nautilus [tT7|, [TBj, would succeed in measuring the impinging cos- 
mic rays in correlation with the resonator mode. Such a result, as a test for the 
further applicability of the thermo-acoustic conversion model at operating temper- 
ature, would come closest to the real situation envisaged for the new detectors. 
Apart from such temperature effects, the applicability of the thermal acoustic con- 
version model [j4|, |^, is confirmed by the data and therefore cosmic rays should be 
expected to seriously disrupt, as calculated by the model, the possibihty of detecting 
gravitational waves. It is beyond the scope of this article to go into any detail [JT^ 



We want to point, however, to earlier calculations 0, |7|, which, having used the 
model, clearly show, firstly, that a next generation spherical resonant mass gravita- 
tional wave detector of ultra high sensitivity will be significantly excited by cosmic 
rays. Secondly, the high impact rate of cosmic rays will prohibit gravitational wave 
detection at the earth's surface, with the required sensitivity. 



Finally, shielding the instrument by an appreciable layer of rock as available in, for 
instance, the Gran Sasso laboratory, would suppress the cosmic ray background by 
a factor of > 10^. Even then a vetoing system will be necessary and, at the radical 
reduction of the background rate so established, it may indeed work effectively. 
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Appendix: Sphere excitation model calculation 



Our calculation of the (L, M)-mode excitation strengths is based on the source term 
ofeq. 5.10/11 of ref. s = J:/{pV)* J dzV±-u, with J: = -fdE/dz. Here, 7 is the 
Grueneisen constant, pV de sphere's mass, and dE/dz the absorbed energy per unit 
track length. The Fourier amplitudes, measuring the second time derivative of the 
mode amplitudes, are directly proportional to s, and the mode amplitudes follow 
from s/cj^, as in eq. 5.18 of [jT^. However, the amount of energy absorbed per unit 
length in our case depends on the particle's position along the track. We therefore re- 
included the E-term under the source term's integral by letting dE{z)/dz represent 
the electromagnetic cascade development of ref. as an approximation to the 
amount of energy absorbed per unit track length by the sphere at position z along 
the beam track, 

f ^ / .dEiz) , , , 

SL,M^ = Kj^y±- UL,MAZ)—^dz, (11) 

where z is being measured from the beam's entrance point into the sphere. With Eabs 
being the total amount of energy absorbed by the sphere from the electron bunch, 
we write dE{z)/dz = Eats * d{E{z) / Eabs)/dz and use the polynomial expansion 
d{E{z)/Eabs)/dz = Yl\^QCiz\ Then J^^"^^'^ dzd{E{z) / Eats) /dz=l. For the polynomial, 
measuring z in meter, we acquired the values co=0.8332 m^^, ci=226 m^^, C2=- 
1832 m~^, C3=4909 m~^ from a fit to the form given in ref. with less than a 

percent deviation for our case of < z < 0.15 m. The value for the energy absorbed 
by the sphere from a single electron, £"^^^=123 MeV, we got from both our Monte 
Carlo simulation using GEANT , and from EGS4 |23, |2|] . At the measured 25 nC 



beam pulse charge this corresponds to a total Eabs = 3.1 J absorbed by the sphere. 
Then the value of k Eabs = jEabs/M=1.00 m^/s^, for our case of M = pV=4:.95 
kg and 7=1.6. Our sphere has a suspension hole which leads to a slight shift in 
the frequencies and the spatial distribution of the modes, with respect to those of a 
sphere without a hole []25[. We approximated, however, our sphere's modes by the 
ideal hole- free sphere's eigenmode solutions u{z) using the available computer 
code as established in ref. [)27j], and renormalising to / u-udV = V, as used in ref. [0 
from eq. 5.6 onward. The source term sl,Ml was calculated for each mode {L,M) 
by numerically integrating eq. [TT]. We checked the surface term in the numerical 
procedure to be negligible, as assumed in the partial integration leading to the form 
of s used by ref. [|T^. Each u(0) in eq. [T^ is the eigenmode solution, calculated 
for each sensor on the 0-grid of the measured data, and each term sl^MlM 
excitation factor sl^Ml the specific beam position bk- 



